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To make certain quantitative interpretations of spectra from NMR experiments carried out on 
heterogeneous samples, such as cells and tissues, we must be able to estimate the magnetic and 
electric fields experienced by the resonant nuclei of atoms in the sample. Here, we analyze the 
relationships between these fields and the fields obtained by solving the Maxwell equations that 
describe the bulk properties of the materials present. This analysis separates the contribution 
to these fields of the molecule in which the atom in question is bonded, the "host" fields, from 
the contribution of all the other molecules in the system, the "external" fields. We discuss the 
circumstances under which the latter can be found by determining the macroscopic fields in the 
sample and then removing the averaged contribution of the host molecule. We demonstrate that 
the results produced by the, so-called, "sphere of Lorentz" construction are of general validity in 
both static and time- varying cases. This analytic construct, however, is not "mystical" and its 
justification rests not on any sphericity in the system but on the local uniformity and isotropy, i.e., 
spherical symmetry, of the medium when averaged over random microscopic configurations. This 
local averaging is precisely that which defines the equations that describe the macroscopic fields. 
Hence, the external microscopic fields, in a suitably averaged sense, can be estimated from the 
macroscopic fields. We then discuss the calculation of the external fields and that of the resonant 
nucleus in NMR experiments. 



INTRODUCTION 

Overview 

NMR spectroscopy is notable for its contributions 
to the study of the chemical and physical properties 
of heterogeneous samples including living cells and tis- 
sues. Variations in the magnetic characteristics of a 
sample often bring about readily observable changes 
in resonance frequency and spectral line shapes, thus 
providing unique probes of cellular function (e.g., 1- 
4) . A knowledge of the physics of systems with mul- 
tiple compartments of differing magnetic susceptibil- 
ity has already laid the foundation for new sorts of 
NMR experiments. The insightful article by Chu et 
al. (5) explains some fundamental aspects of contrast 
enhancement in magnetic resonance imaging (MRI) 
that are brought about by paramagnetic metal-ligand 
complexes; and, the comprehensive review in this jour- 
nal by Levitt (6) gives an independent explanation of 
some key phenomena. Both articles emphasize the 
nature and value of the magnetic field "experienced" 
by a nucleus in a magnetically polarizablc medium. 
They use the theoretical construct of the, so-called, 
"sphere of Lorentz" in their analysis. However, this 
theory provoked us into some deeper questions that 



*Email: p.kuchel@mmb.usyd.edu.au 



seemed to warrant exploration; and, with this insight 
came the expectation of a better understanding of ex- 
perimental data and new experimental methods. 

The accompanying (preceding) article illustrates 
the results from some simple practical NMR exper- 
iments in which samples were chosen in which there 
were differences in magnetic susceptibility across the 
(micro) boundaries of the heterogeneous samples. The 
experiments were conducted on solutions, an emulsion 
in the presence of a solution of the same substance, 
and a suspensions of red blood cells (RBCs) made rel- 
atively paramagnetic. 



Fields in NMR 

NMR spectroscopy is based on the interaction of 
the spin and magnetic moment (s) of a nucleus with 
the magnetic field in its neighborhood. To simulate 
NMR experiments, the magnetic field in the immedi- 
ate vicinity of the nucleus in the "host" atom must 
be calculated. The atom itself may be free or bonded 
as part of a molecule or in rapid exchange between 
these states. In what follows, we refer to the nucleus 
as residing in a host molecule with the understand- 
ing that "molecule" should be interpreted as simply 
the host atom when it is not chemically bonded. An 
NMR spectrometer has a magnet that is designed to 
produce a strong uniform field, which we shall refer to 
as the applied field, into which the sample is placed. 
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As a result of its introduction, the field in and around 
the sample is perturbed by the interaction of the field 
with the magnetic moments of the molecules in the 
sample. Most molecules possess no intrinsic magnetic 
moment because the electrostatic binding forces lead 
to a net cancellation of the orbital and spin angu- 
lar momenta of the electrons. In the presence of an 
applied field, the orbital moments precess and gener- 
ate an extra component that is aligned opposite to 
the field; this is the diamagnetic effect. On the other 
hand, a few atoms, ions, and molecules, in which the 
angular momenta of the electrons do not cancel, pos- 
sess an intrinsic dipole that is much larger than the in- 
duced moment at room temperature. These molecules 
tend to adopt the least-energy state in which the in- 
trinsic dipole moment is parallel to the applied field; 
this is the paramagnetic effect. The intrinsic paramag- 
netic dipole is much larger than the induced diamag- 
netic dipole at room temperatures, so paramagnetic 
molecules are often introduced as "agents" to make a 
deliberate modification to the applied field (e.g., 5). 
However, the magnetic field experienced by a nucleus 
in the host molecule is modified by the fields produced 
by all the molecules that are external to the host in 
which it is located and by the field produced by the 
host molecule itself. The former field is sometimes 
called the local field but because it is in the immedi- 
ate external environment of the host molecule we shall 
refer to it as the external field. The field from the host 
molecule we shall refer to as the host field. The exter- 
nal field is dependent on the composition and geom- 
etry of the whole sample and is a macroscopic entity, 
whereas the host field depends solely on the structure 
of the particular molecule and is a microscopic entity. 
The change in resonance frequency due to the former 
is called the bulk magnetic susceptibility shift (BMS) 
and that due to the latter is the chemical shift. (See 
our preceding article for illustrations of chemical and 
BMS shifts that occur in samples in containers with 
spherical and cylindrical geometries.) 

Both effects arc the result of the electromagnetic 
properties of charges in motion and are therefore de- 
scribed by the Maxwell equations (e.g., 7, 8). In mi- 
croscopic form, these equations describe exactly the 
electric field e and the magnetic induction b pro- 
duced by each constituent-charged elementary par- 
ticle. Hence, these equations can be used to calcu- 
late the host field and the chemical shift at a nu- 
cleus. However, it is not realistic to solve the array 
of Maxwell equations for all the molecules in the sam- 
ple, let alone all those in the microscopic system in the 
vicinity of the nucleus of interest. Therefore, it is nec- 
essary to invoke a simplification that uses the Maxwell 
equations for the macroscopic system, and then we 
work down to the microscopic system. The macro- 




FIG. 1: Figure 1 Schematic representation of a solution of 
pentagonal molecules with nuclei denoted by open discs. 
The relative orientation of the molecules is random and the 
nucleus of interest is placed at the arbitrarily chosen ori- 
gin of a Cartesian coordinate system. The three regions 
of space that are relevant to the averaging, or smooth- 
ing, process used in analyzing the fields experienced by 
the nucleus at the origin are denoted host, external, and 
macroscopic. 



scopic equations describe approximately the bulk elec- 
tric field E and the bulk magnetic induction B, given 
the constitutive relationships between the electric field 
and the electric displacement D, and between the mag- 
netic induction and the magnetic field H. These rela- 
tionships incorporate the effects due to the paramag- 
netic and diamagnetic properties of the bulk medium, 
described above. But, the macroscopic field calcu- 
lated using the macroscopic Maxwell equations does 
not provide the average field at a point where the nu- 
cleus resides within the host molecule. The macro- 
scopic field at any point in fact contains an average 
contribution from the host molecule itself based on the 
assumption that the point is distant from the host 
molecule, not within it. This average host molecule 
field contribution must be removed from the macro- 
scopic field to estimate the average field across the 
molecule due to the applied field and the contribu- 
tions of all the other molecules, the external (or lo- 
cal) field. The actual field experienced by the nucleus 
is obtained by adding the external field component 
and the internal contribution from the host molecule, 
namely, the microscopic host field (see Fig. I). 
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Analytic Strategy 

The commonly used prescription for estimating the 
field at the nucleus of interest (the host field) is based 
on the idea of the sphere of Lorentz (e.g., 5, 6, 9). It 
was introduced in electrostatics by Lorentz (10) but 
applies equally well to magnetostatic situations (e.g., 
11). It was first used in the context of NMR by Dick- 
inson (12). The sphere of Lorentz is a notional sphere 
drawn around a nucleus; it is large enough for all the 
molecules external to it to be treated as a macroscopic 
continuum that is locally uniform. Within the sphere 
the host molecule is imagined to reside in a vacuum 
surrounded by individual molecules whose net elec- 
tromagnetic effect is taken to be vanishingly small. 
The effect of introducing a spherical cavity into a uni- 
form continuous medium is significant and yet the field 
arising from it can be readily calculated; hence, this 
provides the estimate of the external field. 

The disadvantage of this approach is that it mixes 
the microscopic and macroscopic pictures and the jus- 
tification for the field estimates at each level is not en- 
tirely clear. It has thus assumed an almost "magical" 
air; for example, Springer (9) writes "For an imagi- 
nary object, the sphere of Lorentz produces amazingly 
profound real effects. 



Aims, Approaches, and Outcomes 

In this article, we set out to demystify the esti- 
mation of the external fields from the macroscopic 
field(s). We show that the results obtained by us- 
ing the sphere of Lorentz argument are in agreement 
with our new more rigorous approach and are of gen- 
eral validity. This is important because the sphere of 
Lorentz argument is simple to apply and can yield a 
semiquantitative description of BMS effects in vari- 
ously shaped objects. 

Our new analysis requires only that the measure- 
ment "polls" a sufficiently large number of molecules 
that are in random positions with respect to their 
neighbors for the sample average to be equivalent to 
the average of a single molecule surrounded by a ran- 
domized medium. Then, local isotropy, i.e., spherical 
symmetry, produces the same result as the sphere of 
Lorentz construction for a nucleus residing at its cen- 
ter in the otherwise empty space. 

We first summarize the derivation of the expressions 
that describe the macroscopic fields, which are the 
average of the effect of all molecules present in the 
sample; this averaging takes place over macroscopic 
length and/or time scales. On the microscopic scale, 
all charged particles are taken to be in motion so that 
the electric and magnetic properties are coupled, so 



we must analyze the full set of Maxwell equations. It 
is only on the macroscopic scale that experiments can 
realize time scales that are sufficiently long that the 
macroscopic equations governing the magnetic field 
may be solved independently of those governing the 
electric field. At this stage we confine our discussion to 
the magnetic field in the slowly time- varying situation. 
Next, we examine how the average fields produced 
by all the other molecules at the site of a particular 
molecule can be estimated by using the macroscopic 
equations. The method of solution of the macroscopic 
equations is summarized in the appendix; hence, the 
calculation of the BMS shift in NMR experiments is 
completed. 

The detailed treatment of the chemical shift pro- 
duced by the internal field of the host molecule at any 
of its nuclei is beyond the intended scope of this arti- 
cle (e.g., 9, 12, 13) but, for completeness, we provide 
an estimate of the host field in the spirit of the treat- 
ment of the BMS shift so that the relative magnitudes 
of the combined effects can be seen. 

MICROSCOPIC AND MACROSCOPIC FIELDS 
Maxwell Equations 

The Maxwell equations are treated in all standard 
texts on clcctromagnetism (e.g., 7, 8), so we present 
them with little discussion. 

The microscopic electric and magnetic fields e and 
b produced by the moving charges of the submolecular 
particles are described by the following equations: 

<9b 

e V ■ e = p, V x e = - — 

at 

V-b = 0, —V xb = j + fo| [1] 
Mo oi 

where the constants eo and po are the electrical per- 
mittivity and magnetic permeability of free space, re- 
spectively. We must treat the microscopic electric and 
magnetic fields together in a coupled system because 
the microscopic charge density p and current density 
j vary on short time scales due to the rapid motion of 
the point-like charges. 

The Maxwell equations are written here in their dif- 
ferential form using vector differential operators. The 
properties of these operators and their application in 
electrostatics arc described in (14). But, their integral 
form is more readily visualized. Two arc obtained by 
integrating the two equations on the left over a volume 
V and transforming the left sides to integrals over the 
surface S of V, using Gauss's theorem (12): 

J J (e • n)dS = i JJJ pdV, Jj ' (b • n)dS = [2] 
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FIG. 2: Illustrations of Maxwell's laws. Left: The electric 
flux is defined as the integral of the component of the 
electric field e in the direction of the normal n to the 
surface S. It is proportional to the total charge within the 
volume V enclosed by the surface S. Right: The magnetic 
flux crossing any surface S enclosing a volume V always 
vanishes. 



These are illustrated in Fig. 2. The unit vector nor- 
mal to the surface S at any point is denoted n so 
that the scalar product, b • n, is the component of 
b in the direction normal to the surface. The first 
equation thus states that the normal component of e 
integrated over the surface, the electric flux across the 
surface, is proportional to the total charge contained 
within the volume V. This is Gauss's law. The second 
equation states that the total magnetic flux across the 
surface surrounding any volume V must vanish. The 
inward flux exactly balances the outward flux. This is 
due to the empirical fact that there are no sources of 
magnetism corresponding to point charges. 

The other two forms are obtained by integrating the 
two equations on the right in Eq. [1] over a surface S 
that is bounded by a closed curve C, and then trans- 
forming the left side to line integrals around C, using 
Stokes's theorem (14): 

je-d\ = -j t JJ (b-nJdS, 
±-jh-d\ = n)dS + e j t H(e- n)dS [3] 

Here, dl represents a small vectorial increment in the 
path around C. These are illustrated in Fig. 3. The 
first equation states that the line integral of the elec- 
tric field e around a closed path C, known as the elec- 
tromotive force, is equal and opposite to the rate of 
change of the magnetic flux across the surface S. This 
is Faraday's law of induction. The second equation 
states that the line integral of the magnetic field b 
around a closed path C is governed by two quanti- 
ties. The first term on the right side is the current 
crossing the surface S, i.e., threading C. This term 
expresses Ampere's law. The second term is propor- 
tional to the rate of change of the electric flux across 
S. This term was introduced by Maxwell, who called 



FIG. 3: Illustrations of Maxwell's laws. Left: The elec- 
tromotive force is the integral of the component of the 
electric field e in the direction of the vector path incre- 
ment dl around the curve C. This integral is equal and 
opposite to the rate of change of magnetic flux across the 
surface S. Right: The corresponding integral of the mag- 
netic field b around C is determined by the rate of change 
of the electric flux across the surface S and the total cur- 
rent threading the surface. 



it the displacement current, and is essential for the 
existence of electromagnetic waves. 

These equations can be solved by eliminating the 
electric field to derive the equation governing the mag- 
netic field and vice versa. The result in each case is 
a wave equation with source terms provided by the 
electric charge and current densities: 



V 2 < 



V 2 b- 



1^1 

c 2 dt 2 
1 d 2 b 



-- (-Vp- 
-/i V x j 



c 2 dt 



[4] 
[5] 



where c = 1/^/eoMo is the speed of light. The wave 
equation can be solved formally in terms of an inte- 
gral, over all space and all time, of the sources on the 
right sides of these equations. They admit both "ad- 
vanced" and "retarded" solutions, i.e., one in which 
the effect of the source propagates backward in time 
with the speed of light and one in which it propagates 
forward in time. The advanced solution is discarded 
to avoid violating causality. The resulting retarded 
solution can then be expressed as an integral over all 
space of the contribution from every source at the re- 
tarded time. This time is simply the time at which 
the electromagnetic signal, traveling at the speed of 
light, must have left the source point to arrive at the 
selected point of space x where the fields are measured 
at the selected time t. 

This formal integral solution can be written in a 
variety of different ways. The Jefimenko form (e.g., 
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7) is one of the most illuminating: 



Temporal and Spatial Averages 



e(x,t) = ^/// (X : XV ^'W 



47re 



+ 



47T60 /// 



dy(*',t')/dt' 



Cr X — X' 



|x - X'| J 
(x-x')(dp'(x',i')M') 



C X — X' 



l\2 



d 3 x' 



[6] 



Mx,*) = * fff j(x 'f )x( ;: x ° A' 

47T ./././ X-X' 



t /|3 

Mo fff (()j'(x',i')/9i')x(x-x') A , [?] 



Without the final term, Eq. [6] is Coulomb's law and, 
again without the final term, Eq. [7] is the Biot- 
Savart law. In these equations, t = t — |x — x'|/c is the 
retarded time referred to above, and the integration 
is taken over all of space. In the steady state, only 
the first terms of the right side of Eqs. [6] and [7] 
appear, so that the electric field is seen to depend 
only on the charge distribution, and the magnetic field 
depends only on the current distribution. The size of 
the second term relative to the first is l:(d/cr), where 
d and r are the characteristic distance and time scales 
of the system. If the system changes only slowly (viz., 
over time scales long compared with d/c, the time for 
light to travel a distance d), the second terms and 
the variation in the retarded time can be neglected; 
we call this situation the quasisteady case. At the 
microscopic level, the time scale is determined mainly 
by the orbital speeds of the electrons that are an order 
of magnitude less than c. At the macroscopic level 
in NMR experiments, the time scale is governed by 
diffusion, i.e., the thermal speeds of the molecules and 
these are many orders of magnitude smaller than c. 
The quasisteady approximation is, therefore, good in 
the latter case, but more care needs to be taken at the 
microscopic level. 

We now focus on the conceptual steps involved in 
proceeding from the microscopic picture in which the 
molecular constituents of the material are described as 
charged particles in motion in a vacuum, for which the 
solutions can be written down exactly, to the macro- 
scopic picture, which involves systematic averaging 
over the sample. Only by doing this is the relationship 
between the two pictures revealed. Because there are 
several steps that are mathematically independent but 
logically connected, we use a notation that, although 
somewhat cumbersome, does help keep track of the 
steps. 



The development presented by Jackson (7) is based 
on Russakoff (15) and we follow his approach here: 
The Maxwell equations are linear so they can be aver- 
aged over space or time to describe the average micro- 
scopic electric and magnetic fields e and b . These are 
expressed in terms of the sources of the fields, which 
are the charge density p and current density j that are 
averaged in the same manner, namely, 



e V • e 



Vxe + - = 



1 de - 

Vb = 0, — Vxb-eo— =j 
Mo at 



[8] 



The form of these equations is precisely that of Eq. [1] 
so that the integral forms and the general solution ap- 
plied to the averaged fields in the sources are replaced 
by the averaged sources. 



Significance of Linearity 

The property of linearity enables the fields to be ob- 
tained by averaging the sample point over a volume in 
space or interval in time; this is precisely the value of 
the field that would be obtained if it were measured 
at a fixed point but arose from sources that were aver- 
aged over an equal volume or equal time interval. As 
a further consequence of linearity, the average fields 
and sources of any system are equal to the sums of 
the average fields §j and due to the average sources 
p and j representing each of the N molecules present 
in the system. Specifically, e = Thus, each 

average molecular field is described by the Maxwell 
equations, Eq. [8], with the corresponding average 
molecular sources. 

Suppose that the ith molecule consists of rij charged 
particles, namely, electrons and protons; it also con- 
tains neutrons that generate no electromagnetic ef- 
fects per se. We shall not be interested in following 
the submolccular motion, so we first average over a 
time scale that is long compared with that of the or- 
bital motion of the electrons but short compared with 
that of the motion of some reference point (usually 
taken to be the center of mass) of the molecule. We 
denote any time-averaged value of a property Xi with 
an overbar: 



Xi(-K,t) = J Xi(x,t - T)w(r)dT 

where the temporal smoothing function w(t) is nor- 
malized so that J w(r)dT = 1 and it vanishes for time 
scales longer than that of the submolecular motion. 
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i-th molecule 



x 



FIG. 4: Vectors used to specify a general position in 3- 
space, x, the position of the center of the ith molecule, 
Xj(£), and the position of the jth electron in the ith atom, 
£ij(t), relative to the molecule's center. 



By treating each submolecular particle as a discrete 
point charge, the charge and current density due to 
the jth charge in the molecule can be expressed in 
terms of a delta function (which has, it should be 
noted, units of volume -1 ) as 

Pij(x) = qijS(x-Xi(t) - €ij(t) 

& (x) - qa (±i (t) + iij (i)) J(x - x, (t) - & (t)) [9] 

Here, Xj(i) and Xj(t) are the position and velocity of 
the center of mass of the ith molecule; qij is the charge 
of the jth submolecular particle in this ith molecule; 
and £ij(t) and &j(t) are the position and velocity of 
the jth submolecular particle with respect to its center 
of mass (see Fig. 4). 

The delta function <5(x — x') is a convenient mathe- 
matical device for representing a point particle. This 
function vanishes everywhere except at the point 
x = x', yet its integral over any volume contain- 
ing the point x' is unity. Hence, the quantity p = 
q6(x— x') is a charge density such that the total charge 
J I Iv = if the point charge at x' is outside the 
volume V and the total charge J J J v pdV = q if the 
point charge at x' lies within V. Likewise, the quan- 
tity j = <?v<5(x — x') is a current density such that 
the current J J J v jdV = if the moving charge at x' 
lies outside V and the current / / j v jdV = qv if the 
moving charge at x' lies within V . Here, v is just the 
rate at which the point charge is moving; because its 
position vector is a function of time, x' = x'(t), its 
velocity is v = dx'(t)/dt = x'(t). 



Multipole Expansion 

At a point external to the molecule the inequality 
|x— Xj| > holds, so we can expand the expressions 
in Eq. [9] as series in terms of increasing powers of the 
small quantity |£ij|/|x — Xj|. This procedure leads to 
a multipole expansion. Normally, higher-order terms 
are small enough to be negligible compared with the 
first one or two terms. Truncation of the series after 
the second term is referred to as the dipole approxi- 
mation and it yields the following expressions for the 
total charge and current densities for the ith molecule: 

Pi{x,t) = qi S(x-Xi(t)) - V- (pi(^(x- Xi (i))[10] 
d 

ji(x,i) = —(pi{t)5(x - Xi(t)) 

+ Vx(m,(()i(x-x,(i))) 
+ qi±i(t)8(x - Xi(i)) 

+ V X (pi(t) X Xj(t)J(x - Xj(t))) [11] 

where the total molecular charge qi (independent of 
time) , electric dipole moment pi , and magnetic dipole 
moment in, are defined by 

q l = ^2 qij, Pi(*)=5Z 3 «^'W 

= ^E^'^'W x 4-(*)) [12] 

The significance of each of the remaining terms will be 
explained below but, in the meantime, these expres- 
sions suffice to convey the idea that when averaged 
over time the whole molecule appears at distant points 
as if it were a point object. However, it is conceptually 
convenient in what follows (although not mathemat- 
ically imperative) to restore the finite extent of the 
molecule when we deal with its immediate neighbor- 
hood. Thus, we smooth the ith molecule over a volume 
Vi that is comparable to the volume that it occupies. 
The spatially smoothed quantities are denoted by a 
double overbar: 

Xi(x,t) = JJJ x t (x-C,t)h t (C)d 3 C 

where hi describes the form of the smoothing function 
around x. It vanishes outside Vi and is normalized so 
that its integral over Vi is unity. 

The Maxwell equations written with these averaged 
sources describe fields that vary smoothly over length 
and time scales greater than those of the submolecular 
structure. However, the sources and fields still have 
features on the scale of intermolecular distances. The 
macroscopic equations are based on eliminating this 
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fine structure by being smoothed over a volume V of 
space that is sufficiently large to contain a large num- 
ber of molecules. These newly smoothed quantities 
are denoted by angular brackets: 



III 



X t ( x -(,t)h v (()d 3 C 



where hy is another normalized function, like a Heav- 
iside step function, that vanishes outside V. 

The macroscopic averaging process warrants some 
discussion. It is implicit in this concept that the result 
will describe a physical system that can be measured 
with macroscopic equipment and produce repeatable 
results, at least within an acceptable error range. As 
a result, the macroscopic average must be made over 
a volume large enough to ensure that the movement of 
the molecular constituents both within this volume as 
well as in and out of this volume results in little change 
to the average properties. Further, the average should 
reflect the results of measurement by different instru- 
ments that may sample the volume in a similar, but 
not identical, manner. So, the macroscopic aver- age 
should not weight heavily any localized region within 
itself. It is therefore logical, as well as expedient be- 
cause it simplifies the mathematical treatment, to take 
the smoothing function to be uniform over V so that 
hy = 1/V. Thus, we restrict the integration to the 
volume V around x, writing 

(l 4 (x, *)) = 1 HI X 4 (x - C, t)d 3 C [13] 

Combining all smoothings, the average charge and 
current densities become 



(Pi(x,i)) = qiH v (-x.--Xi(t)) 

- V-(pi(t)H v (x-Xi(t)) 

(ji(x,t)> = jt(p.(t)H v (x-Xi(t)) 



[14] 



+ Vx(mi(t)H v (x-Xi(t))) 

+ g i x i (t)fr v (x-x i (t)) 

+ Vx(pi(t) xxi(t)H v (x-Xi(t))) [15] 



where 



ffy(x) = l jjj hitx-QfC 

Thus, the charge sources are contributed by the molec- 
ular charge density and spatial variations of the elec- 
tric dipole density; the current sources are contributed 
by the molecular charge flux, temporal variations of 
the electric dipole density, and spatial variations of 
the magnetic dipole density and the electric dipole 
flux density. 



Finally, the macroscopic charge and current sources 
are obtained from Eqs. [14] and [15] by summing over 
all molecules. Hence, 

(p(x,i)) = -V-P(x,i) [16] 

(j(x,t)) = ^!VxM(x,i)| 

V x ^Pi(i) x ii{t)H v {x-Xi(t))^ [17] 

where the macroscopic polarization P and magnetiza- 
tion M are defined to be 

N 

P(x,t) = ^2pi(t)H v (x-Xi(t)), 

i=l 
N 

M(x,t) = ^mi(t)H v (x-Xi(t)) [18] 

i=l 

At this stage of the analysis, we made the usual as- 
sumptions that the molecular charges and charge 
fluxes qiki sum to zero over a macroscopic volume, 
so the expressions for both the average charge den- 
sity and average current density depend solely on the 
dipole terms. The second term in the expression for 
the average current density is the spatial average of 
the (vector) product of the dipole moment and the 
center-of-mass velocity of the molecules. If there is no 
correlation between the dipole moment and the veloc- 
ity (one is microscopic, the other is macroscopic) then 
the average of the product is equal to the product of 
the averages. In this case, if there is no bulk motion in 
the system the average velocity of the center of mass 
vanishes and the macroscopic current density reduces 
to 



<j(x,t)> = 



0P(x,f) 

Ft 



+ V x M(x,i) 



[19] 



This yields the Maxwell equations that describe, 
within the approximations detailed above, the macro- 
scopic fields E and B in their standard form in the 
absence of free charges: 

<9B 

e V-E = -VP, V x E = - — 

at 

1 dY* <9E 

V • B = 0, —V xB = — + VxM + e — [20] 
fi at at 

These equations should be contrasted with the orig- 
inal microscopic form in Eq. [1]. They have exactly 
the same structure but the source terms are no longer 
discontinuous functions of the microscopic charges and 
of the currents produced as they move. Instead, the 
sources are continuous functions of the electric and 
magnetic dipole densities, which are the highest-order 



8 



terms to survive the averaging processes. However, 
the general solution of these equations is the same as 
for the microscopic equations given by Eqs. [6] and 
[7] with the appropriate change in the expressions for 
the source terms. In the quasisteady case these are 



E(x,t) = - 



1 



47re n 



(x-x')V-P(x',i) - , 



f /|3 



dV[21] 



B(x,t) 



Mo f//VxM(x',/)x(x-x') 



47T 



x — x 
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+ 



In practice, these are just formal solutions because the 
dipole moment densities are, in general, not prescribed 
functions of space and time. Following Maxwell, we 
must therefore proceed by introducing the new fields 



D = e E + P. H = — — M 

Mo 



[23] 



Then, the new Maxwell equations assume the stan- 
dard macroscopic forms: 



<9B 

V • B = 0, — = -V x E 
at 

V • D = 0, ?P = VxH 

at 



[24] 



We now have four vector field quantities, i.e., 12 
scalar quantities, to determine from only eight compo- 
nent equations. This is impossible without additional 
information about how the four fields are related to 
one another. But, before considering this question, 
we can use the prescription for obtaining the macro- 
scopic equations to answer the fundamental question 
of how such fields can be used to estimate the fields 
experienced at the site of a constituent molecule, i.e., 
to calculate what we have defined to be the external 
fields. 



EXTERNAL FIELDS 
Key Concept 

The key concept in understanding the theoretical 
construction derived in this article is: The external 
fields at the site of an individual molecule within a 
sample differ from the macroscopic fields as calcu- 
lated above because the macroscopic field already con- 
tains an averaged contribution to the fields from the 
molecule itself. It is obvious that the macroscopically 
averaged fields experienced by the kth molecule are 
found simply by solving the macroscopic equations, 
with sources from the averaged contribution of that 



molecule subtracted. It is not obvious, however, that 
the actual external fields will always be such fields 
obtained by macroscopic spatial averaging. Indeed, 
in crystals, they may never be. 



Crystals 

In the case of a crystal with molecules arranged in 
a regular lattice, there will be no spatial smoothing 
and the external fields will be similar for all similar 
molecules and could differ greatly from the macro- 
scopic average. The external fields can only be esti- 
mated by solving the full set of microscopic equations 
for the whole lattice; however, this problem is not dis- 
cussed further here as it is not pertinent to molecules 
in solution. 



Amorphous Solids 

In an amorphous solid, with random structure, an 
individual molecule also does not experience spatially 
smoothed external fields. Consider a group of like 
molecules with random arrangements of their neigh- 
bors in a volume V. If this volume is used to define 
the macroscopic average then the external field expe- 
rienced by the group will be similar to the macroscop- 
ically averaged field. The mean field experienced by 
a smaller group (in a smaller volume) will, in general, 
differ from the macroscopically averaged field because 
the system will retain structure on the scale of inter- 
molecular distances. The evaluation of the external 
fields in the latter case poses an intractable problem. 



Averaging Process for Fields and Sources 

The question of how to estimate the mean exter- 
nal field thus hinges on the nature of the averaging 
process. Recall that linearity implies that averaging 
fields is equivalent to averaging the sources. In a fluid, 
the molecular sources move around, so a temporal av- 
erage is equivalent to averaging over the locations of 
the molecules, i.e., to spatial averaging around a fixed 
site. This is the ergodic hypothesis of Boltzmann and 
is plausible, although difficult to prove rigorously in 
most cases. If a finite volume containing a large num- 
ber of molecules with host atoms is sampled, each such 
molecule will be surrounded by other molecules in a 
series of random realizations. The external field, av- 
eraged over such an ensemble of realizations, is again 
equivalent to a spatial average around a fixed site. 
The distribution of surrounding sources is then locally 
uniform and isotropic within the averaging volume V. 



This property characterizes the macroscopic average 
and defines how large the necessary macroscopic av- 
eraging volume should be. We shall consider this case 
exclusively. 

Focus on this averaged spatial distribution of 
sources in the neighborhood of an individual molecule, 
labeled k. To maintain the identity of the molecule, 
we assume that it occupies a volume Vk about the 
"center" of the molecule at x^; all other molecules 
are excluded from this volume. We now smooth the 
molecules outside this volume over a spherical volume 

V to produce a continuous distribution that maintains 
the volume Vk free of sources. As before, we shall take 
the smoothing function to be uniform over V . If V 
does not contain Vk, the sources will be averaged over 
the whole of V so that hy = l/V at points inside 

V and hy = outside. Then, the spatially averaged 
sources contributing to Eqs. [14] and [15] will take the 
form 



N 



X(x,t) = - Jjj^Ux-Ct^t 
X(x,t) = l||^fjl l (x-C,t)d 3 C [25] 

The prime on the sum indicates that the term i = 
k must be omitted. The second expression follows 
from the first because vanishes throughout V. It 
is precisely the definition of the macroscopic average 
(X) in Eq. [13]. 



Smoothing Over Sources in the Neighborhood of 
a Host Molecule 



If we smooth (average) over a volume that includes 
any part of V4 we must confine the sources to the 
volume V minus V£, which is that part of Vk within 
V (see Fig. 5). This volume depends on the position 
of the center, x, of the smoothing sphere relative to 
the center of the molecule x^. The normalization now 
gives hy = 1/(V — VI) and the average sources are 



because X k vanishes throughout V — Vj.. This does 
not match the definition of the macroscopic average 




FIG. 5: Molecular volumes Vk, centered at x fe , and the 
smoothing volumes V, centered at x. Left: The molecule 
lies wholly outside V. Right: The molecule lies wholly 
inside. Center: That part of the molecular volume lying 
inside V is denoted V{.. 



but it can be further expressed as 



X k (x-(,t)d 3 ( 



Now, the first integral is the total source Qy within 
V. Because V is the macroscopic volume, this source 
is given by Eq. [25] as Qy = V(x). The molecules 
will be randomly distributed over this volume with an 
average number density 1/V m , where V m is the aver- 
age molecular "volume" so we can define the average 
source associated with each molecule to be 

Qm = ~^~Qv = Vm{X) 

unless the total source vanishes, i.e., Qy = V(X) = 
0. The second integral is the source due to the fcth 
molecule within V^. By definition, the source within 
V k is 

Q k = JJJ v l fc (x-C,t)d 3 C 
If Xk is uniform within Vk, then it follows that 



X k (x-(,t)d 3 (=^Qk 



If X is not uniform, we can always write, 

J fj^ l fe (x - C, t)d 3 C = /x(x - x fc )^Q fc 



where fx is some smooth function of position with re- 
spect to the center of the molecule and is a function of 
the structure or shape of the kth molecule. When all 



10 



the molecules are of the same type Q k = Q m ', other- 
wise, we write Q k — a k Q m , where a k is a proportion- 
ality constant that depends on the other molecular 
structure(s). Then, substitution yields 



X(x,i) = 



V lit t V * 
— {X - a k f x 



iv-v k ) 



[v-vi)v k 



V! = 

-±(X) 



if (X) ^ 0. This case describes the contributions of 
the molecular dipoles p k and to the sources in 
Eqs. [14] and [15]. When (X) = 0, 



X(x,t) = -f x 



n 



(V Vl)V k 



This case describes the contribution of the possible 
molecular charge q k to Eq. [14] because 

Qk^qk / / / h k (x - x; - ()d 3 C = qk 
J J Jv k 



At locations within V k we must, of course, set 
X(x,t) = 0. 



Alternative Representation of the Sources 

The results of the previous section can be put in an 
alternative, perhaps more illuminating, form by not- 
ing that they differ from the sources for the macro- 
scopic fields only within V k or when V includes some 
part of V k ■ The external fields are therefore obtained 
from the macroscopic fields E and B by removing the 
fields due to the contribution of the fcth molecule to 
the sources within V k , and around V k . We call these 
the self-fields. They are produced by the sources, as 
follows: 



or, if (X) = 0, 



{(X) if x is within V k 

yx(x — x k )(X) if V about x contains any part of V k 
if V about x is wholly outside V k 



{0 if x is within V k 

S^qIx — Xfc)Qfe if V about x contains any part of V k 
if V about x is wholly outside V k 



r 



The shape functions, 

y x = 



OtkfxV m 



(v-v>) 



Vk 



[26] 



contain all the effects of the particular structure of the 
kth molecule. The other quantities are all macroscopic 
quantities. 

The expressions for the self-fields then result from 
solving the Maxwell equations with charge densities: 



{— V • P if x is within V k 

^QQk — V • (J?^P) if V about x contains any part of V k 

if V about x is wholly outside V k 

and current densities: 

r dP/dt + V x M if x is within V k 
j(x, t) = < d(y p P) jdt + V x (y M M) + y Q q k x + V x (,y p P x x) if V about x contains any part of V k 

to if V about x is wholly outside V k 



r 



Prescription for External Fields 



vidual molecule due to a surrounding randomized con- 



The previous section provides a prescription for 
evaluating the external field experienced by an indi- 
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FIG. 6: Successive smoothing processes. Left: The host 
molecule is surrounded by other molecules in random 
positions and orientations. Middle: The surrounding 
molecules have been smoothed into a continuous source 
distribution outside the host molecule. Right: The orien- 
tation of the host molecule has been averaged, resulting in 
a spherically symmetrical distribution of sources outside a 
spherical "exclusion" volume. 



the charge densities are given by 

(- -V • P if r < R k 

p(x, t) = \ y Q {r)q k - V • y p (r)-p if R k < r < R + R k 
I if R + R k < r 

and the current densities are given by 
j(x,t) = 

r dP/dt + VxM if r < R k 

\ d(,y p (r)P)/dt + Vx(,y M (r)M) if R k <r<R + R k 

{ if R + R k < r 

where P and M are the macroscopic polarization and 
magnetization around x = x k . 



Isolated Molecules 



figuration of other molecules. The equations may be 
solved if the structure and dynamics (position, orien- 
tation, and velocity) of the molecule are known. In 
practice, the dynamic state of individual molecules 
will not in general be known, so this formalism has 
limited utility. However, in NMR experiments the sig- 
nal is generated by a large number of nuclei residing 
in molecules, all in different dynamic states. There- 
fore, we can perform another averaging: This time 
it is done over the velocities and orientations of the 
host molecule. If there is no net flux of the molecules 
in question, i.e., there is negligible bulk motion, the 
velocity-dependent terms will average to zero. If we 
also average over all orientations then the newly aver- 
aged sources must be distributed with spherical sym- 
metry. The host molecule is thereby replaced by one 
that is spherically symmetrical, as in Fig. 6. The 
molecular volume V k will then be a sphere and all 
spatially dependent factors in the shape functions y 
will depend only on distance r = |x — Xfc| from the 
center at x k . 

If we recall that the point Xfc is some reference point 
in the molecule, not necessarily the center of mass, it 
is clear that we can now choose that point to be the 
location of the nucleus in the host molecule. When 
the orientations of the molecule are averaged about 
this point, the averaged sources due to all the other 
molecules become spherically symmetrical about the 
nucleus. The average fields experienced by that nu- 
cleus are those due to the spherically symmetrical 
source distributions, evaluated at the center, i.e., at 
x = x fe . 

We can now construct the sources explicitly. Let- 
ting R be the radius of V, and R k the radius of 14, 



Consider the idealized case in which there is a sin- 
gle molecular species (so a k = 1 and q k = 0) and set 
Vk = V m and fx = so that each molecule occupies 
the same exclusive spherical volume that is equal to 
the molecular volume. Then, the shape functions van- 
ish, leaving only the uniform sources P and M within 
V k . Substitution into the quasisteady expressions Eqs. 
[21] and [22] yields explicit estimates of the sclfficlds 
at x = Xfe that are given by 

E. elf = --Lp(x fc ,t) ) B self =^°M(x fe ,i) [27] 

The contribution from the polarization P to B solf in- 
tegrates to zero at the center of the sphere. This does 
not occur at other points within the sphere. How- 
ever, in the strictly time-independent case the fields 
are uniform within a sphere with uniform polarization 
and magnetization so Eq. [27] gives the estimate of the 
static self-fields at all points within the sphere. This 
case reproduces the result obtained from the sphere 
of Lorentz construction, even though the spheres are 
differently conceived and are of different size. The 
contradictions inherent in this simplistic view of the 
sphere of Lorentz are apparent, however, so we pursue 
the rigorous treatment. 

Rigorous Treatment 

If we wish to identify values averaged over V in the 
general case with the macroscopic averages, we must 
take V 3> Vfc. The macroscopic sources are then uni- 
form throughout Vfc and its surroundings, within the 
much larger volume V about the molecule. Thus, in 
the surrounding shell only the shape functions vary 
and then only with respect to the radial coordinate 
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r. The importance of this property was noted by 
(6); substitution of these forms in Eqs. [21] and [22] 
produces exactly the same estimates as those in Eq. 
[27] for the fields at the center of Vk. Away from 
the center, the time-varying polarization will again 
contribute to Bsclf. In the strictly time-independent 
case, the fields will again be uniform throughout Vk, 
so these estimates will apply at all points within Vk- 
The external fields experienced at any position x 
within an averaged molecule in averaged surroundings 
are the macroscopic fields less the self-fields. In the 
static case, these are 



E cxt (x)=E(x) + — P(x), 
B ext (x) = B(x) - ^M(x) 



[28] 



In the time-varying case, these expressions are exact 
at the center of the symmeterized molecule and ap- 
proximate the fields close to the center. 

The most significant conclusion, however, follows 
from the fact that the nucleus of the host molecule 
is placed, by construction, at the center of the sym- 
meterized distributions. The expressions in Eq. [28] 
therefore give the average fields experienced by the nu- 
cleus as a result of the surrounding molecules exactly, 
even if the fields are not static. 



Alternative Expression for the External Field 

An alternative form for the expression for external 
fields can be obtained by recognizing that our assump- 
tions allow us to rewrite Eq. [18] as 



P = Np m , M = Nfh % 



[29] 



where TV is the number density of molecules and p m 
and m are the mean molecular electric and magnetic 
dipole moments in the sample. Then, 



E„ 



„ N _ 

E+ — Pm, 

3e 



B, 



B 



-rci r , 



[30] 



In the case of diamagnctic and paramagnetic 
molecules, which are of considerable interest in NMR 
experiments, the external fields experienced by a 
molecule determine its microscopic (molecular) elec- 
tric and magnetic dipole moments. The mean mo- 
ments of all the molecules in the sample are linearly 
related to the external field by the expressions 



m, 



7m tj 



[31] 



where 7e is the average molecular polarizability and 
7m is the average molecular magnetizability; these 



quantities are provided by the analysis of molecular 
dynamics and have the dimensions of volume. (Note 
that these 7 e , 7 m parameters are not the magnetogyric 
ratio that is usually denoted by this symbol in NMR 
theory) Substituting for the external fields from Eq. 
[30] it is seen that the moments are aligned with the 
macroscopic field: 

_ eo7eE _ ( 7 /Mo)B 

Pm ~l-7eA73' m ~l + 2 7 mA/3 [d2J 

The polarization and magnetization are therefore 

e o7e 7VE 



P = Np m = 
M = Nm m = 



l-7eJV/3' 
(7^/Mo)B 
1 + 2 7m iV/3 



and 



E 



B 



Eoxt l- 7e iV/3' ~° xt l + 2 7m A/3' 



[33] 



[34] 



Susceptibilities 

When calculating macroscopic fields, it is more 
usual to introduce the susceptibilities Xe and Xm, 
which are defined by 



Xe 



7miV 



l-feN/3' 



l- 7m iV/3' 



so that 



D = eE, H = -B 



[35] 



where e = eo(l + Xe) and /i = /io(l + Xm) are the 
permittivity and permeability of the material, respec- 
tively. 

Thus, in terms of the susceptibilities, P and M are 
given by 



P = Xe£oE, M 



Xr. 



B 



(1 + Xm) MO 



hence, 

E ext =(l + f)E, B^=(l- 2 -- 



(1+Xm) 



[36] 



B [37] 



The external fields can now be evaluated directly 
from the macroscopic fields using Eq. [37], where the 
macroscopic fields are found by solving the Maxwell 
equations, Eq. [24], together with the constitutive 
relations, Eq. [35]. 
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Review 

So far, we have developed the theory of external 
fields for both the electric and magnetic fields for sev- 
eral reasons. First, it is based on a model of mi- 
croscopic charges for which the electric and magnetic 
fields are strongly coupled because of the rapid mo- 
tion of the submolecular particles. Second, the macro- 
scopic fields remain coupled, but much more weakly 
because the macroscopic motions are much slower 
than the microscopic ones and produce changes in 
the macroscopic properties only over much longer time 
scales. Third, the external fields can be analyzed sys- 
tematically when both the magnetic and electric fields 
are varying slowly, the analysis leading to the remark- 
able generalization of the sphere of Lorentz construct 
to cases in which time variations are present. The re- 
sults of the construction may therefore be used when 
an electric field is imposed on a sample in an NMR 
experiment or when the applied magnetic field fluctu- 
ates. 

We can therefore safely adopt the simplification 
made in most NMR applications that the system is 
in a macroscopically stationary state because any real 
deviations from this state will not affect the manner 
in which the external fields can be estimated. The ad- 
vantage of the stationary assumption is that the equa- 
tions governing the electric and magnetic fields are 
then decoupled, allowing one to be treated indepen- 
dently of the other. Therefore, in the examples that 
follow to illustrate the theory attention is restricted 
to the static magnetic field. Then, the calculation of 
the external magnetic field from Eq. [37] requires a 
knowledge of the expression for the macroscopic field 
B. This is found via the macroscopic Maxwell equa- 
tions; the time-independent forms are 

V • B = 0, VxH = 0, 

B = M H = Mo (l + X )H = Mo (H + M) 

Their formal solution is treated in standard texts such 
as (7, 8) and an outline in the context of NMR is given 
elsewhere (e.g., 16) and in the appendix. 

HOST FIELD 
General 

The macroscopic fields, and the fields from the 
charged particles, are calculated using the approxi- 
mation in Eqs. [10] and [11], which is appropriate for 
points at a large distance from the system of charges 
that constitute each molecule. When calculating the 
fields experienced by the nucleus of a host molecule, 



there are contributions — the host fields — from the sys- 
tem of charges constituting that molecule. This sit- 
uation requires us to evaluate the fields at an inter- 
nal point of the host molecule. Due to the proxim- 
ity of the electrical and magnetic field sources, these 
fields will in general be more intense than those gener- 
ated by the other, distant, molecules. Their effect on 
the NMR resonance frequency of a nucleus, called the 
chemical shift, is likely to be greater than that of the 
external field, which produces the bulk susceptibility 
shift. The host fields must therefore be determined 
with care for each atomic nucleus in each molecular 
species and it requires the full panoply of quantum 
mechanics. This has been applied by many authors 
from Ramsey (17) onward and is beyond the intended 
scope of this article. However, the effect may usefully 
be illustrated by a simple model. 

Simple Model 

In the absence of an applied field, there will be no 
preferred direction for a molecule and so the averaged 
sources will have a distribution with spherical symme- 
try about any fixed reference point, the nucleus in the 
host molecule, for example. Suppose first that the nu- 
cleus is at the center of a host atom so we assume that 
the positive charge resides at the center. The nega- 
tive charge due to the averaged electron cloud will 
appear as a spherical shell around it. In the presence 
of an external field, the first-order perturbation of the 
spherical shell will be a dipole term with its axis in the 
direction of the applied field. Take the preferred direc- 
tion to be that of the unit vector k. We can take the 
nucleus as the origin of the coordinate system without 
loss of generality. Then, the mean charge density due 
to the electrons about the nucleus can be written as 

Pk (x, t) = po (r, t) - p x (r, t)k • - 

r 

where r = |x| is the distance of the point x from the 
center. The charge density is here a negative quan- 
tity and the center of the distribution is shifted in the 
negative k direction due to the Lorentz force due to 
the local applied field. The mean current density can 
likewise be written 

jfe(x,i) =ji(r,t) (k x -) 

Note that the spherically symmetrical current term 
vanishes and the first-order term is a circular current 
system about the axis k. When j\ is positive this 
represents a paramagnetic effect; when negative, it 
represents a diamagnetic effect. 

The quasisteady microscopic host fields at the origin 
are found by substituting these forms into the appro- 
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priatc terms of Eqs. [6] and [7] and performing the 
integrations, giving 

1 



e fc (t) - ( ^ / pi(r,t)dr)k, 



2/* 



b fc (t) - [~y J ji(r,t)drjk 

We can also evaluate the average electric and magnetic 
dipole moment of these distributions from 

Pfc(t) = - (y/ PiM> 3 dr^k, 
m fc (t) = (y J ji(r,tydr^jk 

The electric dipole is always oriented antiparallel to 
the applied field and the magnetic dipole is parallel 
to the applied field for paramagnetic molecules and 
antiparallel for diamagnetic molecules. Hence, 



efe(i) = 
b fc (t) = 



1 J Pidr 
Air €q J r^pidr 
Mo Jjidr 



Pk(*), 



27T J r 3 jidr 



m k (t) 



The internal electric field is always antiparallel to the 
electric dipole moment and the internal magnetic field 
is always parallel to the magnetic dipole moment. The 
exact relationship between the dipole moments and 
the internal field, of course, depends on the electronic 
structure of the molecule. The result can be written 
in terms of the effective volumes of the molecule V e 
and V m , defined by 

V e = J-p^, V m = 2j-p^ [38] 
J P\dr J ]idr 

e fc (t) = ^rPk(t), b k (t) = -g-m k (t) [39] 

eo V e V m 

Similar expressions will therefore arise in the more 
complicated cases that occur when the host atom is 
chemically bound to a molecule. Averaging all possi- 
ble orientations of the molecule will produce charge 
and current distributions with spherical symmetry 
about the reference point, which is the nucleus of the 
host atom. The external field will produce first-order 
perturbations of these distributions that give rise to 
net dipole moments and to the related host field at 
the nucleus (center) of the host molecule. 



DISCUSSION 

Sphere of Lorentz Argument 

We have shown above that the sphere of Lorentz 
construct does indeed provide a means of estimat- 



ing the fields, both electrical and magnetic, at points 
within a molecule embedded in a macroscopic sample 
composed of other molecules. Our derivation does not 
depend on an ad hoc "hard" spherical construction 
but it can be viewed as sort of "soft" spherical one. 
The analysis demands that we do not attempt to esti- 
mate the fields experienced by an individual molecule, 
but we estimate the average fields experienced by a 
large collection of similar molecules so that both the 
dynamic properties of the molecule and the locations 
of the neighboring molecules are randomized. Under 
these circumstances, the microscopic electromagnetic 
sources in the close vicinity of the molecule display 
spherical symmetry, as do the macroscopic ones. It 
is the assumption of spherical symmetry that pro- 
duces the general results described by Eq. [28], not 
sphericity. This explains why the radius of the as- 
sumed sphere of Lorentz plays no role in the result; 
the properties within a sphere of any radius drawn 
around the center of a molecule will be symmetrical 
and lead to the same result. Our argument is, how- 
ever, not independent of the distance scale so the size 
of the averaging sphere is not totally irrelevant. The 
macroscopic properties of the material are required 
to be locally uniform, so the size of the macroscopic 
smoothing volume must be chosen to ensure this ap- 
parent homogeneity. This volume must be sufficiently 
large to randomize the contribution of the least abun- 
dant molecular species. Therefore, the smoothing vol- 
ume will be least when all the molecules in the sample 
are of the same species. For a macroscopic model to 
make sense, this volume must be less than that of any 
macroscopic heterogeneities. Such cases will be illus- 
trated in the examples below. 



Chemical Shift 

The magnetic field experienced by a nucleus in 
an atom in a molecule is dependent upon the local 
bonded structure of the atom. This field is the ba- 
sis of the experimentally important chemical shift of 
the intrinsic resonance frequency away from the value 
that pertains to the isolated atom. In what follows 
is an explanation of this phenomenon and a model of 
a molecular system that provides an estimate of the 
order of magnitude of the chemical shift effect relative 
to that of the BMS shift. 

It is immaterial whether the host atom, in which the 
nucleus of interest is located, is chemically bonded into 
a larger molecule or not. In cither case, the nucleus 
will experience, on average, fields that are the sum of 
the external fields, Eq. [28], and the host fields, Eq. 
[39] . The size of the host fields will, of course, depend 
on the molecular environment of the host atom. Nev- 
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ertheless, the dipole moments of the atom or molecule 
to which it is bonded are linearly related to the ex- 
ternal fields as in Eq. [31] but with the mean molecu- 
lar polarizability and magnetizability replaced by the 
quantities specific to the type of molecule in question, 
k say. Carrying out these substitutions gives the host 
fields at the nucleus, Eq. [39]: 



efe(i) = - 



le,k 



E 



b fc (t) = 



e 'm 

Then, we find the total fields at the nucleus to be 



E, = E...-e, = (l + ^jh-^* |E 



Bfe B oxt + b fe 



2X, 



3(1 + Xm) 



1 + M| B 



considerations simply emphasize the need to have ac- 
curate estimates of the host field at the site of the 
nucleus of the host atom when the molecule is aver- 
aged over all orientations. In this situation the re- 
quired analysis of chemical shifts is sophisticated and 
computer intensive (e.g., 9). 



ILLUSTRATIVE EXAMPLE 

In the following we illustrate the theory by evaluat- 
ing the external field and the susceptibility shift in a 
suspension of RBCs following the procedure outlined 
by Wolber et al. (18). 



Susceptibility-Induced Shifts 



where E and B are the macroscopic fields. These can 
be put in a more revealing form by writing 7 eJ t in 
terms of the mean polarizability 7 e as 



le.k = 



le,k 



le 7V(l+ Xe /3) 

and similarly for the magnetizability. Then, 



p _ | 1 , Xe V.. 7e,fc 1 Xe 
' ' 3 A 7e W e (l+Xe/3) 



1 



2Xr 



X 1 + 



3(1 +Xn 
Jm,k 1 



Ira NV m (1 + Xm/3) 



B 



[40] 



E [41] 



[42] 



Now, N is the number of molecules per unit volume 
so 1/NV is the ratio of the average intermolecular 
volume to the volume of the molecule, as defined by 
Eq. [38]. In fluids this ratio is, say, of the order of 
10, falling to unity when the molecules are densely 
packed. The other factor is the ratio of the specific to 
the mean value of 7 e or j m . If there is only a single 
species present, this ratio is unity so the chemical shift 
due to the host field contribution will be an order of 
magnitude greater than the bulk susceptibility shift 
due to the external contribution. If the molecules are 
diamagnetic then j m < in the host field and this 
reduces the field experienced by the nucleus; so, the 
effect is referred to as diamagnetic shielding. If the 
molecule is paramagnetic, "/ m > and the nuclear 
field is enhanced. When there are several species of 
molecules present, the host contribution depends on 
the properties of the molecule in question relative to 
the average. It will be greatest when the molecule to 
which the host atom is bonded is strongly paramag- 
netic while the mean is weakly diamagnetic. These 



These authors begin their analysis with a system in 
which a uniform strong magnetic field BO is created in 
a material of susceptibility xo- There is then uniform 
magnetization in the material given by Eq. [36]. If a 
sample with susceptibility Xs is now introduced then 
the new field B will be the sum of the original field 
and the field B due to the change in magnetization 



M' 



XsB 



XoB 



Mo(l + Xs) Mo(l + Xo) 



inside the sample and 



M , = Xo(B_Bo) 

Mi + xo) 



[43] 



[44] 



outside the sample. But, B' also satisfies the following 
equation (see appendix, Eq. [75]) for the case of a 
single surface S: 



B' = Mo M' - £° 

47T 



(AM'(x')-n)(x-x') j2 



- v'|3 



dV [45] 



where S is the surface of the sample. As noted in 
the appendix, these equations do not provide a full 
general description of the field. For a full solution, 
it is usual to solve the Laplace equation, Eq. [77], 
with appropriate boundary conditions, for the scalar 
potential (e.g., 16). However, Eq. [45] does provide a 
useful approximate form of the expression when the 
magnetic susceptibilities are small. In practice, x 1S 
of the order of -10 x 10~ 7 cgs-emu (-40tt x 10~ 10 SI 
units) so that |B'|/|B | is also of order —10 x 10~ 7 . 
Hence, to first order in small quantities we can write 



M' « (Xs Xo) B 
Mo 



[46] 
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inside the sample and M' = outside. Again, to first 
order, 



B' « (x. 



+ 



Xo)B 

(Xs - Xo) 
An 



(B .n)(x-xQ 

I v — -v' I 3 



[47] 



The new macroscopic field inside the sample is thus, 
approximately, 



B(x) = B +(xs-Xo) 



x B 



The integral over the surface of the sample is a func- 
tion solely of the geometry of the surface. In general, 
it is a function of the position x of the point within the 
sample. For ellipsoidal surfaces (including the special 
case of a spherical surface), however, the result is in- 
dependent of x; in other words, the field within the 
sample is uniform. For a sphere, 



^ Us 



(B -n)(x-x') , 2 , B n 



AitJJg |x-x'| 3 
so the macroscopic field is 



B 



l+gfe-Xo )B 



[49] 



[50] 



In general, we may write 



where 2) s is the geometric demagnetizing factor that 
is described by (18) among others; and, for a sphere 
2> s = 2/3. For an infinite cylinder aligned at right 
angles to the field, 2> s — 1/2. Of course, there is no 
surface effect for a cylinder aligned parallel to the field 
and 2> s — 1. A table of such factors is provided by 
Chu et al. (5). The geometric factor is the same at 
all points within these objects and in all objects with 
ellipsoidal surfaces (16); the field is therefore uniform 
within them. For other objects, the geometric factor 
will be a function of position and the field will not be 
uniform. In either case, we can write the macroscopic 
field as 



B = (1 + (x« — Xo)^)B 
and the external field as 

2Xs 



[52] 



B, 



1 



3(1 + Xs) 



B 



Bo + (Xs - Xo)^ s B - -XsB [53] 



to first order in the susceptibilities. This expression 
describes the external field in absolute terms, i.e., it is 
the field experienced at a point surrounded by a local 
vacuum in the sample. The expression agrees with (5) 
but it agrees with Eq. [1] in (18) only if the original 
material is taken to be a vacuum, i.e., xo = 0. 

The difference arises from the focus of the latter 
authors on the frequency shift observed in NMR ex- 
periments. This shift is governed by the change in 
the external field as a result of introducing the sam- 
ple. Before the introduction the macroscopic field is 
B so the external field at the site of a nucleus is ap- 
proximately 



*L = (l--Xo)Bo 



while the external field Bf xt at the site after the in- 
troduction of the sample is given by Eq. [53]. The 
change is therefore 



AB, 



which is Eq. [1] in (18). 



Xo)B 



Sample Heterogeneity 

Wolber et al. (18) also consider the case in which 
the sample is heterogeneous, consisting of plasma and 
RBC. If the magnetic field in the sample is averaged 
over a volume large enough to contain many erythro- 
cytes, the system can be considered to be uniform with 
a volume-average susceptibility given by 



Xb = 



VeXe + VpXp 
V e + V p 



[54] 



where the subscripts e and p denote erythrocyte and 
plasma, respectively, and V e is the average volume of 
an erythrocyte and V p the average volume of plasma 
surrounding each erythrocyte. The macroscopic field 
that determines the external field at a nucleus can 
then be calculated as outlined above from the aver- 
ages bulk properties. If sample size is not large com- 
pared to the size of the heterogeneities in the system, 
account has to be taken of the detailed distribution 
of erythrocytes in the neighborhood of the nucleus in 
question. To do this, these authors conceptually sur- 
round the site by a sphere, £ , that is comparable in 
size to the averaging volume (see Fig. 7). The field 
source outside this sphere contributes at the center 
an approximate field, called the "far" field, that is 
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FIG. 7: Representation of a heterogeneous system of RBCs 
as randomly dispersed spheroids; in reality, mammalian 
RBCs are biconcave discs. The sample is bounded by the 
surface S. In the second panel, the system is treated as 
continuous outside the virtual sphere E of Wolber et al. 
(18). Within the individual erythrocytes have surfaces Si. 
In the third panel, equivalent to the second, the system 
is treated as continuous outside the selected RBC with 
surface So- 



described by 

>/ _ (Xb-Xo) ff ( B o -n)(x-x') 2 , 



B' 



+ - 



4tt JJs |x-x'| 3 
{Xb-Xo) ff (B -n)(x-x') 2 , 



A-k 



X — X 



'13 



dV [55] 



the normal being directed out of the sample in the first 
integral and into the sphere in the second. As a result, 
the two terms will almost cancel one another if the 
averaging volume is comparable to the sample volume; 
however, the purpose of the construction is to place as 
much as possible of the sample outside the surface E 
so that it can be treated in an average manner. Using 
the results from above, these expressions yield 



B' « {Xb - Xs)(@ s - 1)B + (Xs - Xo)^B 



= (Xb-Xo) [ %~ 3 J B 



[56] 



This is the generalization of Eq. [3] of (18) for xo 7^ 0. 
To this must be added the "near" field contributed 
by sources within the sphere; this field is estimated 
by taking an alternative macroscopic average over a 
volume smaller than that of an erythrocyte but larger 
than that of any internal structures. This construction 
ensures that the susceptibilities Xb and Xp are uniform 
within the erythrocyte and the surrounding plasma, 
respectively. The macroscopic system is thus divided 
into a hierarchy of macroscopic systems, each of which 
contribute to the macroscopic field at the site of a 
nucleus. Then, the changes of magnetization within 
the spherical volume, which are given by 



M" 



(Xe - Xo)B m „ _ (x P - Xo)B ^ 



produce an approximate near macroscopic field at a 
point in the plasma of 

B£« (x P -Xo)B 

, (x P -X e ). E £ (Bo^)(x-^) ^ [68] 



47T 



1 ff (xs-Xo)(B -n)(x-x') , 2 , 



+ — 

4?r ./ , /s 



,/|3 



dfx' [59] 



Here, the sum is taken of the integrals over the sur- 
faces Si of all the erythrocytes within the sphere with 
the normal at the surface pointing into the ery- 
throcyte; and, the final term is the integral over the 
surface of the surrounding sphere S with outwardly 
directed normal. The susceptibility xs at points on 
this sphere will vary according to whether the point 
lies in plasma or erythrocyte. However, the sphere 
will, by construction, intersect a large number of ery- 
throcytes and so the integral may be approximated by 
replacing by the average bulk susceptibility Xb on the 
surface. In this approximation 



1 

47T 



(Xs - Xo)(B • n)(x - x') rf2 ^, 



x - x' 



' 1 3 



~ -(Xb-Xo^Bo [60] 

which, of course, cancels the contribution of the sur- 
face of the sphere to the field B'. 

The near contribution to the macroscopic field in 
an RBC is, in the same approximation, given by 

K~ (Xe-Xo)Bo 

(Xe-Xp) ff (B • n )(x-x') 2 , 



+ 



47T 



So 



- v'|3 



dfx' [61] 



+ (x^Xp! j^ff (Bo-mOfr-x') ^ [62] 



47T 



Si 



1 ff (XE-Xo)(B -n)(x-x') 2 , 



+ — 

4tt ./ , /s 



X - X' 
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d 2 x' [63] 



Mo 



Mo 



where So is the surface of the selected RBC and the 
sum is taken over all the others (n in total) in the 
sphere, the normal now pointing out of each RBC. 
The first integral is evaluated at an internal point and 
can be treated as before in Eq. [51], replacing s by @ e , 
which is the geometric factor appropriate to the shape 
of the RBC. The integral over the sphere S is treated 
as before but it appears to have been omitted by (18). 
After making this correction, the total macroscopic 
field at a point in an RBC is 

B « (1 + (xe - Xo) + (Xb - Xo)(@s - 1) 

+ (Xe - X P )(@e - 1))B + ( Xe - Xp)B"' [64] 
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where 



To first order in the susceptibilities, the external field 
experienced by a nucleus in an RBC will be 



B oxt w ^1 + (x e - Xo) + (X6 - Xo)(^ s - 1) 

+ (Xe ~ XpX^e - 1) - \x)j B + ( X e ~ Xp)B"' [66] 

This expression corrects and generalizes Eq. [5] in 
(18). These authors introduced a further term to de- 
scribe the contribution of the particles within a sphere 
of Lorentz drawn about a nucleus within the RBC; 
then, they argued that it vanished "because of sym- 
metry." We have shown above that the contribution is 
in fact included in the estimate of the external field if 
the particles are randomized within the "local" macro- 
scopic volume that surrounds the nucleus, so the final 
result is the same. 

Wolber et al. (18) then claim that B'" vanishes 
due to symmetry. However, in general, there will be 
no symmetry in the distribution of the other RBCs 
about the selected cell. An exact evaluation is then 
impossible. We argue that this contribution can be es- 
timated only if the experiment samples enough RBCs 
(either in volume or time) that the configuration ex- 
ternal to any individual RBC is randomized. Then, 
the structured medium can be replaced by a homo- 
geneous medium with a mean susceptibility given by 
Eq. [54]. In this case 

K~ (Xe-Xo)Bo 



+ 



{Xe-Xb) ff (B -n )(x-x') d 2 , 



47T 



. _ „/|3 



+ (X6 - Xo) ff (Bo -n)(x-x') dV 



47T 



,/|3 



= (Xe " X0)B + {X \ Xb \ @e - 1)B 



47T 



- (Xfe-Xo)gBo 



and the external field becomes 



1 + (Xe - Xo) + (X6 - Xo)(@s ~ 1) 



B 



+ (Xe - XbX^e - 1) - gXejBo [67] 

Of course, this is exactly the same result as that which 
would be obtained by ignoring the first spherical con- 
struction and applying Eq. [75] directly to the surfaces 



that define the selected erythrocyte and the sample 
volume. The selected RBC is, on average, completely 
surrounded not by plasma alone but by plasma con- 
taining other RBCs, so the susceptibility in the sur- 
rounding medium is Xb, n °t Xp- 11 we se t Xo = 0, we 
obtain 



B oxt « ( 1 + X b{@. -1) 



+ (Xe-X&)(^e-l)+3Xe)B [68] 



which corrects the result given by (16), 



B oxt — 

( 1 +Xb (®. - fj + (Xe ~ Xp) (®e ~ |))Bq [69] 



CONCLUSIONS 

In conclusion, we can finally consider the valid- 
ity of the assumption of spherical symmetry for a 
molecule averaged over all its orientations: Exact 
spherical symmetry would result if there were no pre- 
ferred directions. However, applied electric or mag- 
netic fields do provide preferred directions and the 
molecules (especially macromolecules) may be dis- 
torted systematically as a result. These departures 
from spherical symmetry will lead to a second-order 
correction being required to the estimate of the exter- 
nal field. But, the first-order effect estimated above is 
of order X ~ 10~ 7 , so second-order effects will be of or- 
der x 2 ~ 10~ 14 so they will be negligible. The expres- 
sions in Eq. [28] therefore provide accurate estimates 
of the averaged external fields. Hence, we have shown 
that the mathematical constructs presented above are 
well defined and allow the bulk susceptibility shift to 
be calculated on a more rigorous basis than hitherto. 

In practice, an NMR experiment samples a large 
number of nuclei in a macroscopic volume over a 
macroscopic time interval, and each will experience a 
fluctuating environment as the neighboring molecules 
move around. As a result, the molecules will not ex- 
perience a single average field but a spread of field 
strengths. This fluctuation will appear as a broaden- 
ing of the resonance line about the mean value that 
can be calculated as described here. In principle, the 
size of the fluctuations can be estimated from the 
width at half height of the NMR spectral line if the 
geometric factors are known accurately (see previous 
accompanying article). 
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APPENDIX 



The steady-state solution of the Maxwell equations 
that describes the magnetic field in a sample is given 
by the first term in Eq. [22]. Recalling the fact that 
the integration is performed over all space, this expres- 
sion can be transformed into two alternative forms by 
integrating by parts and assuming that the surface in- 
tegral vanishes sufficiently far from the sample nuclei. 
These expressions are 



B 



= -/// 



Mo M(x') x (x-x') ^, 

47T |x — X'l 3 



These two forms are equivalent for arbitrary volumes. 

The first form displays the construction of the field 
from the vector potential 



■/// 



/i M(x') x (x-xQ ^, 

47T |x — X'l 3 




FIG. 8: Heterogeneous system with volumes Vi embedded 
in a volume Vo- The surfaces separating volumes Vi and 
Vj are denoted Sij and the normal vector n^- at points of 
the surface Sij is directed from Vi to Vj, where j > i. 



because the Maxwell equation V • B = guarantees 
that we can write B = V x A. 
The second can be rewritten as 



-B M = H 



Mo 



-Ah 



M <*> [71] 



X — X' 



/|3 



which exhibits the construction from the scalar po- 
tential 



-III 



M(x')-(x-x') dV 



47T X — X' 



/|3 



because the Maxwell equation VxH = guarantees 
that we can write H = — V<3?m- 

If M is everywhere differentiable and the integral 
is taken over all space, the scalar potential can be 
rewritten as 



Vx'-M 
47r|x - x' 



;d 3 x' 



[72] 



which is the solution over all space of the Poisson 
equation 



V 2 $ m = V • M 



[73] 



The Poisson equation is just a form of the wave 
equations [4] and [5] when there is no time depen- 
dence. The solutions can therefore be obtained from 
the solutions of the wave equation, Eqs. [6] and [7], 
by performing the trivial integration over time when 



the integrand has no time dependence. This produces 
Eq. [72]. 

If space contains media with discontinuous distribu- 
tions of M, the field may be determined from either of 
the two equivalent forms of Eq. [70] , but the volume 
integrals need to be evaluated with some care. 

Suppose space is divided into n regions with vol- 
umes Vi within which Mj(x) is differentiable. Let Vi 
and Vj have a common surface Sij, which may have 
zero extent (Fig. 8). Then, we can transform the 
vector potential to obtain a contribution from both 
within each volume Vi and from each surface Sa 



+ 



- 3 = 1 



III 



Mj(x'-) x n ij d 2 x > 



X - X' 



VvxM t (x') fe ; 



X - X' 



Here, the magnetization Mj(x'— ) in the surface inte- 
gral is evaluated just inside the volume Vi, the primed 
sum excludes the term with i = j, and n,j is the out- 
ward normal from Vi to Vj on Sij . 

If the magnetization is uniform within Vi, V x M, = 
within Vi and the volume integral vanishes. If the 
medium is uniform so \i is constant within Vi, then 
V x M, = XiV x H within Vi. In the magnetostatic 
case, the volume integral will again vanish. 

If space is composed only of such regions 



n n „ „ 
j=l 7 = 1 J Jb '3 



Mi(x'-) x n. 



■d 2 x' 
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and the terms can be paired because Sij = Sji and 
iiij = -r\ji to give 



A = 

n n 



i=l j=i+l J Ji >i] 



(M»(x'-) - Mj(x'+) ) x n^ , 



X — X' 



where now M J (x'+) is now evaluated just outside the 
volume Vi. The quantity AMy(x') = M 3 -(x'+) - 
Mj(x'— ) is the jump in magnetization as the bound- 
ary Sij is crossed from Vi to Vj (in the direction of 
iiy) at the point x'. 

The field follows immediately: 



B = 



Mo 

4-7T 



n n 



EE 

i=l j=i+l 



(AM tJ (x / ) x n tJ ) x (x-x') ^, 



X — X' 
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Alternatively, we may work from a similar version 
of the scalar potential: 



n / n 



i=l \ j = l J J 



Under the same conditions as before V • M = and 
the volume integrals vanish, leaving 



n n 

~i^E E 



M 



AM^(x')-n^ v 



47T ^ — ' * — ' JJq.. |x— X' 

1=1 j = I+l J Jb '3 



The magnetic field is now given by 
B = /j M 



n n 



SEE II (AMt3( " r "7,,| (x ' x/) rfV [75] 



- v'|3 



i=l j=i+l u J °*J 



where M is the appropriate magnetization at the point 
x. 

In this form it is easier to implement the bound- 
ary conditions when the materials are diamagnetic or 
paramagnetic. In either case, M • n is the normal 
component of the magnetization at the boundary and 



AMy • n = 



(Xj - Xi) 



Mo(l + Xj)(l + Xi) 



B(x') • n 



because Maxwell's equations require the normal com- 
ponent of B to be continuous across a boundary. 



Hence, 



n n 
l+X 47T ^— ' 

JJs.Ai + XiW 



i—1 



(Xj-Xi) ( B ( x ')' n .j)( x - X '),2 / 

s'„(l+Xj)(l+Xi) " |X-X'|3 



dV [76] 



This appears to be an equation that defines B implic- 
itly everywhere but it cannot be implemented as such 
because it does not define the discontinuous tangential 
component of B fields on the surfaces Sij . The gen- 
eral method of solution when V • M = everywhere 
is to solve Eq. [73] with vanishing right side: 



V 2 $ 



M 







[77] 



This is Laplace's equation, a version of Poisson's equa- 
tion without sources. Because there are no sources, 
the construction in Eq. [72] is not applicable and 
the solution of Laplace's equation must be deter- 
mined from the conditions imposed on the bound- 
aries. The boundary conditions on <fr M at the surfaces 
between the different media are fully determined by 
the Maxwell equations, namely, that $m is continu- 
ous across the boundary so that the components of H 
parallel to the boundary are continuous and the com- 
ponent of B normal to the boundary is continuous, 
i.e., (1 + x)( n ' V$m) is continuous. 

These prescriptions allow the magnetic field B to 
be found exactly in all circumstances. In practice, 
however, the smallness of the susceptibilities allows 
approximations to be made that greatly simplify the 
estimation of the field. This is demonstrated in the 
illustrative example in the text above. 
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